Destruction of Majorana modes in a topological nanowire by disorder in the 

superconducting substrate 
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We study the proximity effect between a disordered superconductor substrate and a topological 
nanowire coupled via tunnel coupling. We illustrate how a finite pair-breaking term in the super- 
conductor affects the local density of states in the nanowire. On one hand we show that the tunnel 
amplitude between the superconductor and the nanowire induces a topological transition in which 
the Majorana fermionic states can be destroyed when the coupling is very strong. On the other hand 
we note a non-trivial effect of the finite pair-breaking which can take the semiconducting nanowire 
between a topological and a non-topological phase in certain ranges of parameters. Our results have 
direct consequences for a nanowire coupled to an inhomogenous superconductor. 

PACS numbers: 73.20.-r, 73.63.Nm, 74.45. +c, 74.50.+r, 



I 



i 

C 

o 
o 



CN 

> 

o 

CN 

o 

CO 



X 



Introduction Majorana fermionic states have received 
a lot of interest because of their exotic properties [2, Q 
such as non-Abelian statistics that open the perspective 
of using them for quantum computation 0, 13] ■ A sys- 
tem which is expected to exhibit such states consists of a 
semiconducting wire with strong spin-orbit coupling such 
as InAs or InSb [5( , in the presence of an applied Zeeman 
field and in the proximity of an s-wave superconductor 
(SC) [|1 Q. The properties of this system have been ex- 
tensively studied for an homogeneous sc [Ml. 

A few recent experiments have reported the observa- 
tion of zero-bias peaks floM^ll ] which are in good agree- 
ment but non completely consistent with the existence 
of Majorana fermions. Such experimental observations 
have been shown to also have simpler explanations be- 
cause a zero bias peak could also arise without any Ma- 



jorana fermions 22-26]. Most theoretical studies until 
now have focused on clean superconductors which gen- 
erate a proximity effect in the semiconducting nanowire. 
The effect of an inhomogenous SC on the Majorana states 
has received less attention 
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In this letter we consider a more general approach al- 
lowing one to characterize the proximity of a supercon- 
ductor in the presence of disorder. Thus, we study the 
proximity effect induced by a hopping term between the 
wire and the SC substrate, and analyze the differences 
between this model and the more widely used one in 
which the proximity effect is modeled by an on-site SC 
gap induced in the wire. We apply our model to an s- 
wave superconductor with a finite pair-breaking [3l| . We 
study the effect of the latter as well as of the tunneling 
rate between the superconductor and the semiconducting 
nanowire on the formation of Majorana states. We show 
that increasing the tunnel coupling induces a topological 
transition which splits the zero-energy Majorana peak; 
interestingly enough, this peak can in certain conditions 
reform by increasing the pair-breaking strength, but this 



happens only for a very finely tuned regime of parame- 
ters, and it is unlikely to have observable experimental 
consequences. We also observe that, as expected, with in- 
creasing the pair-breaking strength the Majorana peaks 
widen and get dissolved in the bulk. Our model allows 
us also to characterize the physics of a semiconducting 
nanowire in contact with an inhomogenous superconduc- 
tor. 




FIG. 1. (Color online) Scheme of our system: a semiconduct- 
ing nanowire (Bare Green's function Gq(u>)) connected via 
proximity effect to a dirty superconducting bulk (Supercon- 
ducting self energy Es(w)). 

Model Our starting point is a semiconducting wire with 
strong spin-orbit coupling, in presence of a Zeeman field, 
and in proximity with an s-wave superconductor (see Fig. 
[1]). The total Hamiltonian for such a system is of the form 



H = H 



NW 



Ha + Hn 



(1) 



In terms of the extended Nambu spinors the annihi- 
lation operator of the nanowire electrons is given by 



At - 



— (ci. 



i^p cj , C[., — Of), and the semiconducting nanowire 
Hamiltonian becomes 



Hnw 



i 2 - 



(i)t z + apa y r z + V z a z 



cdx, (2) 



where a and t are the Pauli matrices respectively in spin 
and particle-hole spaces, fj, is the chemical potential, a is 
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the strength of the Rashba spin-orbit coupling and V z is 
the applied Zeeman magnetic field. The Hamiltonian of 
the bulk superconductor can be written as 



(3) 



fc,(7 



,2 

with £fc = j— — and ^fc^ is the annihilation operator 
for an electron in the superconductor, having spin a and 
momentum k. The hopping term between the SC and 
the nanowire takes the following form 



(4) 



where the operator c\ a creates an electron with spin a 
on site i. Since the total Hamiltonian is quadratic in 
the SC degrees of freedom, we can integrate out these 
modes, such that the effect of the SC substrate is taken 
into account by dressing the bare Green's function of 
the nanowire by a superconducting self-energy G^{uj) = 
Gq R (uj) — £fj(w). The total retarded superconducting 
self-energy can be written as Sfj(w) — I® T,f R (u>), with 

X being the unity matrix in the space of sites, and Sf R 
the on-site retarded superconducting self-energy given by 



(5) 



The local self energy depends on the tunneling amplitude 
ti as well as on the retarded Green's function for the 
superconducting bulk electrons and can be written as 



VA2 



(6) 



where 1 is the unity matrix in spin/particle-hole space, 

l ~ l 2 

and Ti : s = Tru(0) ItA is the tunneling rate. We take this 
to be uniform along the nanowire for most of the rest of 
the paper, T it s = Ts, U = t. 

The retarded bare Green's function of the semicon- 
ducting nanowire electrons [34] is given by 



G-^H = (u + iS)l-H t 



NW 



(7) 



where S is an infinitesimal quasiparticle inverse lifetime 
which is introduced to avoid divergences in the numerical 
evaluations. For the purpose of our analysis, the Hamil- 
tonian of the nanowire is best described using a lattice 
tight-binding model, such that Eq. ((2|) becomes 

N 



HNW = E C i t(* ~~ ^ Tz + V z a z} °i 
i 

- \ c \ l tT z + iao y T z + h.c] Cf+i, (8) 
where we defined the creation operator c] of an electron 



( c lf c li' c 4' ~ Ci t)- Here and throughout the remainder of 
the text, we work with physical dimensions corresponding 
to h = 1 , and in units in which the hopping term between 
sites is t = 1. 

This general formalism allows us to model also a disor- 
dered as well as an inhomogenous superconductor. Such 
properties are encoded in the superconducting self en- 
ergy which is related to the SC Green's function. Thus 
we first use our model to study the effect of a finite pair- 
breaking which is taken into account via a pair-breaking 
parameter 7, phenomenologically introduced as an en- 
ergy imaginary part [3~0 ] , ui — > lj + ij in Eq. ([5]) . 

In our calculations we focus on obtaining the local den- 
sity of states (LDOS) (evaluated on a given site i) in the 
wire, which is given by the imaginary part of the dressed 
Green's function 



i( W )=-- E Im[G^(o 

77- < J 



(9) 



/9=t,4- 



in the nanowire on site i in the Nambu basis as c. 



Throughout this letter we focus on the limit where 
the system is in the topological phase by chosing V z = 
0.4, A = 0.3, fi = 0, and a = 0.2. A finite width (5 = 
0.002hvF/a) is introduced in the numerical evaluations 
yielding a finite width of the peaks in the LDOS. We 
consider also that we are in the regime I << £ where £ 
is the superconducting coherence length and / the length 
of the nanowire. 

Induced superconducting gap and the formation of Ma- 
jorana modes We focus first on the qualitative features of 
the proximity effect induced by the tunnel coupling with 
the SC substrate, and on the differences with the model 
in which the pairing term is put "by hand" directly in 
the nanowire 0, E3, EH . As described in Fig. [5| a) , 
in the latter case a superconducting gap A is induced 
in the nanowire (in absence of Zeeman field and Rashba 
spin-orbit coupling). When we turn on the Zeeman field 
the effective gap becomes A — V z , and for V z equal to 
A the gap closes. For values of V z larger than A and in 
the presence of spin-orbit coupling, the gap reopens and 
is topological in nature, and two Majorana states form 
at zero energy. The Rashba spin-orbit coupling localizes 
the two Majorana states more and more at the two ends 
of the nanowire. 

When the proximity effect is modeled via a hopping 
term, in the absence of a Zeeman field and of Rashba 
spin-orbit coupling, a gap-like feature also appears at an 
energy A, which can be described as a SC pseudo-gap 
(see Fig. [5]b) ). However, inside this gap we note the 
formation of Andreev bound states (ABS), the number 
of these states being determined by the ratio between the 
inverse nanowire length l/l and the size of the gap (i.e. 
the energy difference between two ABSs is proportionnal 

to 1/0 [ag. 

In the absence of spin-orbit coupling and magnetic 
field, these ABSs have a uniform weight along the wire 
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FIG. 2. LDOS of two models: DOS dependence on energy 
when the superconductivity is put "by hand" (a) and when 
the superconductivity is induced via a tunnel coupling (b). 



(see Fig. [3] a) ) and are a reminiscence of the quantized 
modes of the wire in the normal state (35j . Turning on the 
spin-orbit coupling localizes these states towards the ends 
of the wire (see Fig. [3] b) ), even in the non-topological 
state. Including a Zeeman field breaks the spin degener- 
acy and splits each of these states. The two ABSs come 
closer and closer with increasing the magnetic field un- 
til when, for a certain value of V Zl and in the presence 
of spin-orbit coupling, they merge (note this merging is 
also controlled by the chemical potential) (see Fig.[3]c) ) 
and form two Majorana states. Increasing the spin-orbit 
coupling localizes these states more and more towards 
the ends of the wire (see Fig. [3] d) ) . 

Thus we see that the proximity effect induced via the 
tunnel coupling between the SC and the nanowire also 
allows for the formation of localized Majorana modes, 
which can be seen evolving from the Andreev bound 
states in the wire. 

Topological transitions In order to understand how 
these Majorana states are formed in a more quantitative 
manner, we analyze the effective Hamiltonian induced in 
the wire via the coupling with the superconducting sub- 
strate. Thus, we note that if we focus on the low-energy 
sector where we expect the Majorana modes to form, 
the self-energy induced in the wire is equivalent with an 
effective SC gap which can be evaluated by setting the 
energy to zero in the second term in Eq. ([6]), yielding 
A c ff = TgA/ a/ A 2 + 7 2 . It is this effective gap which 
governs the transition to a non-topological phase when 
the Zeeman field becomes smaller than a critical value, 
K 2 <A c 2 ff + M 2 . 

Alternatively seen, for large values of A e gp a transition 




FIG. 3. (Color online) ABS dependence on energy and po- 
sition a) non-topological state, no Zeeman field and no spin- 
orbit coupling; b) non-topological state with a finite spin-orbit 
coupling; c) topological state, not too large spin-orbit cou- 
pling; d) topological state, strong spin-orbit coupling. The 
pair-breaking parameter is considered negligible, 7 = 0. 

to a non-topological phase occurs. In terms of the tunnel 
coupling between the SC and the wire the correspond- 
ing condition to exit the topological phase is given by 

r s > T s , where f s = (V z 2 - ^ 2 )(1 + £). This can be 
checked numerically, and in Fig. 2] a) we have plotted the 
LDOS at one end of the nanowire as a function of the en- 
ergy and Ts- Indeed, we note that for values of Ts larger 
than Ts the Majorana peaks split in two, marking the 
exit from the topological phase. We have checked that 
the energy difference 8E between these ABSs depends 
linearly on Tg. 

We also note that for very small tunnel coupling pa- 
rameters the Majorana states are also split; this behavior 
can be explained by the fact that in this regime the super- 
conductivity has not sufficiently penetrated the nanowire 
to have a real superconducting gap. By increasing the 
nanowire length the splitting of the Majorana states oc- 
curs for smaller and smaller values of Ts, confirming that 
this is a finite size effect; for large nanowire lengths we 
obtain a robust zero bias peak even for very small trans- 
mission rates. 

While we do not detail it here, this formalism allows 
one to consider also a non-uniform tunneling rate which 
corresponds to disorder at the interface between the wire 
and the superconducting substrate. Such disorder yields 
a position dependent effective superconducting gap in the 
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wire, and the Majorana modes are destroyed for large 
disorder strengths. 



system to behave as a topological wire in the non-SC 
("normal") regime (no Majorana modes). 
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FIG. 4. (Color online) LDOS of a Majorana peak (on the first 
site of the wire) as a function of energy and the tunneling rate 
Ts (a) and pair-breaking strength (b). The dashed line de- 
notes the transition between a topological and non-topological 
state which occurs at Ts ~ 0.4 (with our parameters). In b) 
we take a Ts = 0.3 in the main plot (no topological transi- 
tion with varying 7). In the inset we take A = 0.15, V z = 1, 
and Ts ~ 1.01 (a Majorana peak is expected to form when 
7 > 7 w 0.021). 

We study also the effect of pair-breaking in the SC 
substrate. In terms of the pair-breaking parameter 7 
the condition to have a topological phase is given by 



7 > 7 = A 



1. Thus, the position of the tran- 



sition depends essentially on the difference T s — V z . For 
large magnetic fields V z > Ts the topological condition 
is satisfied for all values of 7, and no transition to a non- 
topological state is possible. In Fig. |4] b) we have plot- 
ted the LDOS corresponding to a Majorana peak (i.e. 
the LDOS on the first site of the wire) as a function of 
energy and pair-breaking strength for this configuration 
(V z >Ts)- Note that the Majorana peaks widen and get 
dissolved in the bulk when the pair-breaking strength 
becomes of the same order of magnitude as the SC gap. 
This is not very surprising as for such strong pair break- 
ing the BCS DOS is smoothed out and becomes similar 
to that of a normal metal; in this regime we expect our 



When Ts > V z the Majorana fermions are destroyed 
at 7 = 0, however, for certain parameter values chosen 
such that 7 is small with respect to the topological gap 
(Ts very close to Ts, large V z with respect to A) there 
exists a transition to a topological phase with increasing 
7. Finding the range of parameters in which this can be 
realized requires fine tuning, and in the inset in Fig. 4 b) 
we have illustrated such a situation (A = 0.15, V z = 1, 
Ts = 1-01), for which a phase transition exists and the 
Majorana peak reforms for 7 > 7 0.021 

Inhomogenous junction Last, but not least, we focus 
on an inhomogenous superconductor consisting of re- 
gions with different tunnel coupling between the wire and 
the SC, as well as with different pair-breaking strengths. 
In particular we consider a superconductor containing 
three regions, a central region in the topological phase 
(Ts < Ts, 7 = 0), and two exteriors ones for which the 
pair-breaking strength as well as the tunnel coupling can 
be tuned such that the regions are in either the topo- 
logical or the non-topological phase. When all regions 
are in topological phase, we restore a well-known case 
where we only have two Majorana states localized at the 
ends of the nanowire. When only the central region is 
in the topological phase, two Majorana fermions form at 
the extremities of the central part. An interesting sit- 
uation is that of two exteriors regions with very strong 
pair-breaking. In this kind of configuration, the exteriors 
regions are quasi-normal, and we recover the behavior 
described in Ref. 14] for SN junctions, i.e., two Majo- 



rana fermionic states extended over the entire "normal" 
parts, as well as over a broad energy range. 

Conclusions To summarize, we have used a general 
approach to describe the proximity effect induced by a 
bulk superconductor in a semiconducting nanowire via a 
tunnel coupling. We have applied our method to study 
the effect of the disorder in a superconductor (modeled 
as a pair-breaking term) on the physics of the formation 
of the Majorana fermions in the wire. We have shown 
that the nanowire can reach a non-topological phase for 
strong values of the tunneling rate between the SC and 
the wire. The pair-breaking strength can also take the 
system trough a topological phase transition in a finely 
tuned range of parameters. Finally we have shown that 
Majorana states can form in an inhomogenous supercon- 
ductor. Our approach can be also applied to other kinds 
of disorder such as static or magnetic disorder. 
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